ART. 

The  University  of  Texas  at  Austin 


TECHNICAL  REPORT 


13  March  2000  Copy  Number  _ 

ARL-TR-00-01 

ADS  MNT  and  FET: 

Array  Element  Localization  Analysis  and  Results 

Technical  Report  under  Contract  N00039-96-D-0051 ,  TO  No.  0236, 

TD  No.  1406046,  ADS  Acoustic  Test  and  Evaluation  and  System  Engineering  III 


Peter  E.  McCarty,  Marten  I.  Garay,  Mary  R.  Parse,  and  Russell  E.  Henrichs 


20010716  011 


Prepared  for:  Space  and  Naval  Warfare  Systems  Command  Headquarters 
4301  Pacific  Coast  Highway  (OT-1)  •  San  Diego,  CA  92110-3127 

Monitored  by:  Naval  Sea  Systems  Command 
Department  of  the  Navy  ‘Washington,  DC  20362-5101 

Acknowledgment  of  Support  and  Disclaimer  (DFARS  252.235-701 0): 

This  material  is  based  upon  work  supported  by  the  Naval  Sea  Systems  Command 
under  Contract  No.  N00039-96-D-0051 ,  Task  Order  No.  0236. 

Any  opinions,  findings  and  conclusions,  or  recommendations  expressed  in  this  material  are  those  of  the  author(s) 
and  do  not  necessarily  reflect  the  views  of  the  Naval  Sea  Systems  Command. 


Approved  for  public  release;  distribution  is  unlimited. 

Copyright  ©  2000  •  The  University  of  Texas  at  Austin  •  Applied  Research  Laboratories  •  Reproduction  and  Redistribution  Prohibited  Without  Prior  Express  Consent 

Applied  Research  Laboratories^The  University  of  Texas  at  Austin»Post  Office  Box  8029»Austin,  Texas  78713-8029 


UNCLASSIFIED 


REPORT  DOCUMENTATION  PAGE  Approved 

OMB  No.  0704-0188 

Public  reporting  burden  for  this  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searchino  exlstino 
data  sources,  gathenng  and  maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information.  Send  comments  regarding  this  burden  estimate 
or  any  other  aspect  of  this  collection  of  information,  including  suggestions  for  reducing  this  burden,  to  Washington  Headquarters  Services,  Directorate  for  Information 
Operations  and  Reports,  1215  Jefferson  Davis  Highway,  Suite  1204,  Arlington,  VA  22202-4302,  and  to  the  Office  of  Managementand  Budget,  Paperwork  Reduction 
rroject  (U7U4-01 88),  Washington,  DC  20503. 

1.  AGENCY  USE  ONLY  (Leave  blank)  2.  REPORT  DATE 

13-03-2000 

3.  REPORT  TYPE  AND  DATES! 

1 

SOVERED 

echnical 

4.  TITLE  AND  SUBTITLE 

ADS  MNT  and  FET:  Array  Element  Localization  Analysis  and  Results 

Technical  Report  under  Contract  N00039-96-D-0051,  TO  No.  0236,  TD  No.  1406046, 
ADS  Acoustic  Test  and  Evaluation  and  System  Engineering  III 

5.  FUNDING  NUMBERS 

N00039-96-D-0051 

TO  No.  0236 

TD  No.  1406046 

6.  AUTHOR(S)  ~  ~  - - 

McCarty,  Peter  E.  Garay,  Marton  1.  Parse,  Mary  R.  Henrichs,  Russell  E. 

7.  PERFORMING  ORGANIZATION  NAMES(S)  AND  ADDRESS(ES)  - - 

Applied  Research  Laboratories 

The  University  of  Texas  at  Austin 

P.O.  Box  8029 

Austin,  Texas  78713-8029 

8.  PERFORMING  ORGANIZATION 
REPORT  NUMBER 

ARL-TR-00-01 

9.  SPONSORING/MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 

Space  and  Naval  Warfare  Systems  Command  Naval  Sea  Systems  Command 
Headquarters  Department  of  the  Navy 

4301  Pacific  Coast  Highway  (OT-1)  Washington,  DC  20362-5101 

San  Diego,  CA  92110-3127 

10.  SPONSORING/MONITORING 

AGENCY  REPORT  NUMBER 

11.  SUPPLEMENTARY  NOTES 

(See  reverse  side  for  "Acknowledgment  of  Support .  . .  (DFARS  252.235-701  OV  and  for  ARLrUT  convrinht  .<5tflfempnt  i 

12a.  DISTRIBUTION/AVAILABILITY  STATEMENT 

Approved  for  public  release;  distribution  is  unlimited. 

12b.  DISTRIBUTION  CODE 

13.  ABSTRACT  (Maximum  200  words) 

Accurate  array  element  localization  (AEL)  is  an  important  component  of  optimizing  the  detection  performance  of 
bottom-mounted  acoustic  arrays.  Use  of  the  nominal  element  positions  can  lead  to  significant  degradations  in 
beamforming  performance  when  the  array  has  been  subject  to  real-world  deployment  methods.  Methods  of  accurately 
estimating  array  element  positions  using  both  overt  and  covert  methods  have  been  developed  by  the  Advanced 
Deployable  System  (ADS)  program,  and  were  used  during  the  ADS  Fleet  Evaluation  Test  (FET)  and  Multi-Node  Test 
(MNT).  "FET”  henceforth  is  used  to  denote  both  the  FET  and  the  MNT.  During  FET,  data  were  collected  for  AEL 
analysis.  This  report  first  explains  the  AEL  algorithms  used,  and  then  presents  AEL  results  for  the  18  arrays  deployed 
during  FET.  The  algorithm!  generally  achieves  an  accuracy  on  the  order  of  0.5  meters  when  using  acoustic  data  with  a 
bandwidth  of  150  Hz. 

14.  SUBJECT  TERMS 

ADS  /  Advanced  Deployable  System  AEL  /  array  element  localization 

FET  /  Fleet  Evaluation  Test  MNT  /  Multi-Node  Test 

15.  NUMBER  OF  PAGES 

91 

16.  PRICE  CODE 

17.  SECURITY  CLASSIFICATION  18.  SECURITY  CLASSIFICATION 

OF  REPORT  OF  THIS  PAGE 

UNCLASSIFIED  UNCLASSIFIED 

19.  SECURITY  CLASSIFICATION 

OF  ABSTRACT 

UNCLASSIFIED 

20.  LIMITATION  OF 

ABSTRACT 

SAR 

NSN  7540-01-280-5500 


UNCLASSIFIED 


v^iaiiuaiu  rum  I  ££70  \nKSV.  ^-O^) 

Prescribed  by  ANSI  Std.  239-18 
298-102 


11.  SUPPLEMENTARY  NOTES  (cont'd) 

Contract  No.  N00039-96>D-0051 
Acknowledgment  of  Support 

This  material  is  based  upon  work  supported  by  the  Naval  Sea  Systems 
Command  under  Contract  No.  N00039-96-D-0051,  Task  Order  No.  0236. 

Contract  No.  N00039-96-D-0051 
Disclaimer 

Any  opinions,  findings  and  conclusions,  or  recommendations  expressed  in  this 
material  are  those  of  the  author(s)  and  do  not  necessarily  reflect  the  views  of  the 
Naval  Sea  Systems  Command. 

Copyright 

Copyright  ©  2000  •  The  University  of  Texas  at  Austin  •  Applied  Research 
Laboratories  •  Reproduction  and  Redistribution  Prohibited  Without  Prior  Express 
Consent 


11 


TABLE  OF  CONTENTS 


LIST  OF  FIGURES .  v 

PUBLICATION  NOTES .  ix 

1.  INTRODUCTION .  1 

2.  CORRELATION  APPROACH  TO  AEL .  3 

2. 1  TIME  DELAY  MEASUREMENTS  FROM 

OMNIxOMNI  CORRELAGRAMS .  3 

2.2  MODELING  OF  TIME  DELAYS,  AND 

ADJUSTMENT  OF  ELEMENT  POSITIONS .  5 

2.3  ADS  AEL  ALGORITHM  VARIANTS  . .  8 

2.4  QUALITY  ASSESSMENT  OF  ALGORITHM  RESULTS .  14 

2.5  ARRAY  APERTURE  ESTIMATION  USING 

SOURCES  OF  OPPORTUNITY .  17 

3.  IMPLEMENTATION  IN  PAS .  21 

4.  AEL  RESULTS  FOR  FET  AND  MNT .  23 

5.  DISCUSSION  OF  AEL  RESULTS .  33 

5. 1  COMPARISON  OF  RESULTS  AND 

ACCURACY  ESTIMATES .  33 

5.2  AUTOMATION  OF  THE  AEL  PROCESS . 35 

5.3  THREE-DIMENSIONAL  AEL .  37 

APPENDIX  A .  39 

APPENDIX  B .  59 


iii 


This  page  intentionally  left  blank. 


IV 


LIST  OF  FIGURES 


Figure  Page 

2.1  Correlation-based  AEL  algorithmic  flowchart .  4 

2.2  Sample  omnixomni  correlagram  with  track  indicated 

by  red  outline .  6 

2.3  AEL  algorithm  variants .  9 

2.4  Effects  of  applying  far-field  AEL  algorithm  to 

a  mid-field  source  of  opportunity .  12 

2.5  Source-receiver  geometry  upon  successful 

algorithm  convergence .  15 

2.6  Source-receiver  geometry  upon  unsuccessful 

algorithm  convergence .  16 

2.7  Aperture  estimation  for  node  1 1 ,  using 

sources  of  opportunity .  18 

3. 1  AEL  data  flow  in  the  analysis  processor .  22 

4.1  AEL  results  for  MNT-only  nodes .  24 

4.2  AEL  results  for  string  1 .  25 

4.3  AEL  results  for  string  2 .  26 

4.4  AEL  results  for  string  3 .  27 

4.5  AEL  results  for  string  4 .  28 

4.6  Array  aperture,  LOB,  and  curvature  summary .  30 

5. 1  Omnixomni  correlagrams  for  far-field 

and  mid-field  contacts .  34 

B .  1  Node  1 1 .  Omnixomni  correlagram  for  near-field  contact 

during  AEL  event,  with  processed  track  highlighted .  61 

B.2  Node  11.  Navigation  geometry  during  AEL  event,  and 

source  track  upon  convergence .  62 

B.3  Node  12.  Omnixomni  correlagram  for  near-field  contact 

during  AEL  event,  with  processed  track  highlighted .  63 

B.4  Node  12.  Geometry  during  AEL  event  (top)  and 

source  track  upon  convergence  (bottom) .  64 

B.5  Node  21.  Omnixomni  correlagram  for  near-field  contact 

during  AEL  event,  with  processed  track  highlighted .  65 

B.6  Node  21 .  Geometry  during  AEL  event  (top)  and 

source  track  upon  convergence  (bottom) .  66 

B.7  Node  22.  Omnixomni  correlagram  for  near-field  contact 

during  AEL  event,  with  processed  track  highlighted .  67 


V 


Figure  Page 

B.8  Node  22.  Geometry  during  AEL  event  (top)  and 

source  track  upon  convergence  (bottom) .  68 

B.9  Node  23.  Omnixomni  correlagram  for  near-field  contact 

during  AEL  event,  with  processed  track  highlighted .  69 

B.  10  Node  23.  Geometry  during  AEL  event  (top)  and 

source  track  upon  convergence  (bottom) .  70 

B.  1 1  Node  24.  Omnixomni  correlagram  for  near-field  contact 

during  AEL  event,  with  processed  track  highlighted .  71 

B.12  Node  24.  Geometry  during  AEL  event  (top)  and 

source  track  upon  convergence  (bottom) .  72 

B .  1 3  Node  3 1 .  Omnixomni  correlagram  for  near-field  contact 

during  AEL  event,  with  processed  track  highlighted .  73 

B .  1 4  Node  3 1 .  Geometry  during  AEL  event  (top)  and 

source  track  upon  convergence  (bottom) .  74 

B .  1 5  Node  32.  Omnixomni  correlagram  for  near-field  contact 

during  AEL  event,  with  processed  track  highlighted .  75 

B.16  Node  32.  Geometry  during  AEL  event  (top)  and 

source  track  upon  convergence  (bottom) .  76 

B.17  Node  33.  Omnixomni  correlagram  for  near-field  contact 

during  AEL  event,  with  processed  track  highlighted .  77 

B.18  Node  33.  Geometry  during  AEL  event  (top)  and 

source  track  upon  convergence  (bottom) .  78 

B.  19  Node  34.  Omni  correlagram  during  AEL  event, 

with  processed  track  highlighted .  79 

B.20  Node  34.  Geometry  during  AEL  event  (top)  and 

source  track  upon  convergence  (bottom) .  80 

B  .2 1  Node  4 1 .  Omnixomni  correlagram  for  near-field  contact 

during  AEL  event,  with  processed  track  highlighted .  81 

B.22  Node  41.  Geometry  during  AEL  event  (top)  and 

source  track  upon  convergence  (bottom) .  82 

B.23  Node  42.  Omnixomni  correlagram  for  near-field  contact 

during  AEL  event,  with  processed  track  highlighted .  83 

B.24  Node  42.  Navigation  geometry  during  AEL  event, 

where  a  green  "+"  indicates  time  processed,  a  black  "+" 

indicates  time  not  processed,  and  a  red  "o"  indicates 

receiver  location .  84 

B.25  Node  43.  Omnixomni  coffelagram  for  near-field  contact 

during  AEL  event,  with  processed  track  highlighted .  85 


VI 


Figure  Page 

B.26  Node  43.  Geometry  during  AEL  event  (top)  and 

source  track  upon  convergence  (bottom) .  86 

B.27  Node  44.  Omnixomni  correlagram  during  AEL  event, 

with  processed  track  highlighted .  87 

B.28  Node  44.  Geometry  during  AEL  event  (top)  and 

source  track  upon  convergence  (bottom) .  88 


vii 


This  page  intentionally  left  blank. 


PUBLICATION  NOTES 


Acknowledgment  of  Support  and  Disclaimer  (DEARS  252.235.7010) 

Contract  No.  N00039-96-D-0051 
Acknowledgment  of  Support 

This  material  is  based  upon  work  supported  by  the  Naval  Sea  Systems  Command  under 
Contract  No.  N00039-96-D-0051,  Task  Order  No.  0236. 

Contract  No.  N00039-96-D-0051 
Disclaimer 

Any  opinions,  findings  and  conclusions,  or  recommendations  expressed  in  this  material 
are  those  of  the  author(s)  and  do  not  necessarily  reflect  the  views  of  the  Naval  Sea 
Systems  Command. 

Copyright 

Copyright  ©  2000  •  The  University  of  Texas  at  Austin  •  Applied  Research  Laboratories  • 
Reproduction  and  Redistribution  Prohibited  Without  Prior  Express  Consent 


IX 


This  page  intentionally  left  blank. 


X 


1. 


INTRODUCTION 


Between  February  and  May  of  1999,  the  Advanced  Deployable  System  (ADS) 
program  deployed  a  total  of  18  arrays  off  the  southern  California  coast  as  part  of  the 
U.S.  Navy's  Multi-Node  Test  (MNT)  and  Fleet  Evaluation  Test  (FET)  exercises.  Data 
from  these  arrays  arrived  via  underwater  cable  at  a  shore  processing  site,  where  the  data 
were  beamformed  and  displayed  for  detection  analysis. 

Such  analysis  requires  that  the  beamformer  and  postbeamforming  processors  be 
provided  with  array  configuration  information,  including  the  relative  positions  of  each 
receiver  element  within  the  array  and  the  overall  orientation  of  the  array.  Historical 
evidence  indicates  that  deployment  of  a  linear  array,  fully  elongated  and  properly 
positioned,  is  the  exception  rather  than  the  norm  in  array  deployment.  Accurate 
estimation  of  array  configuration  and  orientation  is  necessary  to  ensure  optimal 
performance  of  both  the  individual  nodes  and  the  field  as  a  whole. 

To  ensure  this  level  of  performance,  methods  and  tools  have  been  developed  by 
the  ADS  program  to  perform  array  element  localization  (AEL).  The  correlation-based 
AFT,  algorithm  implemented  at  FET/MNT  evolved  from  methods  developed  over  several 
years  at  ARL:UT,  where  they  were  used  for  analysis  of  many  diverse  data  sets  collected 
by  several  different  programs.  These  methods  were  implemented  and  refined  within  the 
ADS  processing  and  automation  system  (PAS),  primarily  within  the  analysis  processor 
(AP). 
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During  FET  (henceforth  used  to  include  both  the  FET  and  MNT  exercises),  data 
were  collected  for  AEL  analysis.  This  included  data  collected  during  events  involving 
USN  ships  RV7  and  RVl  in  motion  near  each  node  (overt  method)  and  data  collected 
during  events  of  opportunity  (covert  method).  These  data  were  analyzed,  primarily  on¬ 
site  in  quasi-real  time,  but  also  post-test.  This  report  summarizes  the  results  of  these 
analyses. 

In  Section  2,  an  overview  of  the  AEL  process  is  presented,  with  implementation 
details  discussed  in  Section  3.  Results  from  all  events  are  then  shown  in  Section  4,  with 
observations,  accuracy  estimates,  lessons  learned,  and  conclusions  covered  in  Section  5. 
Appendices  list  the  best  available  element  positions  for  each  array,  as  well  as  showing  the 
omnidirectional  (onmixomni)  correlagrams  and  geometries  used  during  analysis. 
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2. 


CORRELATION  APPROACH  TO  AEL 


The  basis  of  the  correlation  approach  to  array  element  localization  is  the 
measurement  of  the  timings  of  a  set  of  signals  at  the  different  receiver  elements.  These 
timings  are  modeled  using  assumed  positions  for  the  receiver  elements.  The  assumed 
element  positions  are  iteratively  adjusted  to  minimize  the  mean  square  error  between  the 
measured  and  modeled  timings,  using  a  Hessian-based  nonlinear  optimization  routine. 
Detailed  discussion  of  this  process  follows. 

2.1  TIME  DELAY  MEASUREMENTS  FROM  OMNIXOMNI 
CORRELAGRAMS 

Traditionally,  timing  measurements  used  for  AEL  analysis  were  obtained  from 
impulsive  events  (explosive  charges  or  lightbulb  implosions),  or  from  matched  filtering 
of  modulated  waveforms  such  as  hyperbolic  frequency  modulations  (HFMs).  Because  of 
the  necessity  for  covert  operation  of  the  ADS  system,  as  well  as  for  the  ease  of 
conducting  exercises  without  any  need  for  source  employment,  an  approach  was 
developed  to  obtain  the  necessary  timing  measurements  directly  from  the  radiated  noise 
of  surface  ships. 

In  this  approach,  illustrated  in  Fig.  2.1,  an  omnixomni  correlagram  is  formed 
between  a  particular  pair  of  receiver  elements.  During  the  transit  of  an  acoustic  source,  a 
trace  will  be  observed  on  the  correlagram.  The  trace  is  caused  by  the  broadband  energy 
naturally  radiated  by  the  ship,  and  will  change  in  bearing  as  the  source  transits.  The  time 
delay  measured  from  the  correlagram  represents  the  difference  in  travel  time  from  the 
source  to  each  receiver  of  the  pair.  The  absolute  travel  time  is  not  required  in  AEL 
analysis. 
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Measure  time  delay  vs  time  along  a  CPA  event  from  an 
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Between  10  and  30  specific  times  are  selected  during  the  transit  of  the  source;  at 
these  points  travel  time  difference  is  measured.  In  addition,  measurements  are  obtained 
for  up  to  3N  omnixomni  correlagrams,  where  N  is  the  number  of  receivers  in  the  array. 
In  this  manner,  a  highly  redundant  set  of  time  delay  measurements  is  obtained,  saved  in 
an  array  of  dimension  [(number  of  source  points)  by  (number  of  omnidirectional  pairs)]. 

Obtaining  this  set  of  measured  time  delays  is  one  of  the  more  difficult  steps  in  the 
AEL  process,  both  in  terms  of  expenditure  of  analyst  time  and  in  terms  of  complexity  of 
automation.  A  partially  automated  program  has  been  developed  to  obtain  time  delay 
measurements  from  a  series  of  correlagrams,  but  the  program  currently  requires  analyst 
supervision,  both  to  correctly  identify  the  trace  of  interest,  and  to  ensure  that  the  tracker 
does  not  jump  to  an  adjacent  trace.  Figure  2.2  shows  an  example  omnixomni 
correlagram  from  FET,  with  the  desired  trace  outlined  in  red.  Other  traces,  some  very 
close  in  bearing,  offer  ample  opportunities  for  introduction  of  errors  by  an  unsupervised 
tracking  program. 

2.2  MODELING  OF  TIME  DELAYS  AND  ADJUSTMENT 

OF  ELEMENT  POSITIONS 

The  next  step  in  the  array  element  localization  process  is  to  model  the  array  of 
time  delays  based  on  assumed  positions  for  the  source  and  receivers.  Using  methods 
discussed  in  Section  2.3,  these  source  and  receiver  positions  are  initialized,  and  then 
travel  time  differences  are  computed,  based  on  this  model. 

Next  comes  computation  of  the  RMS  error  between  these  modeled  travel  time 
differences  and  those  measured  from  the  omnixomni  correlagrams.  An  iterative 
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nonlinear  optimization  algorithm  is  employed  to  adjust  unknown  parameters  until  the 
root  mean  square  (RMS)  error  is  minimized.  In  this  case,  the  unknown  parameters  are 
the  x-y  positions  of  each  receiver  element  (assumed  to  have  constant  depth),  and 
parameters  describing  the  source  positions  or  track.  This  algorithm  iterates  until  the 
RMS  error  decreases  to  less  than  a  preset  tolerance,  at  which  point  the  algorithm  is 
considered  to  have  converged. 

Two  nonlinear  optimization  methods  have  been  employed.  The  first  is  the  inverse 
Hessian  method,  in  which  it  is  assumed  that  the  error  surface  is  quadratic  with  respect  to 
the  unknown  parameters,  and  the  update  to  the  unknown  parameters  is  based  on  a  local 
estimate  to  this  quadratic  surface.  This  assumption  is  generally  good  when  unknown 
parameters  are  close  to  their  optimal  values.  In  terms  of  the  AEL  problem,  this 
assumption  is  good  when  the  error  has  a  simple  dependence  on  the  unknown  parameters, 
and  when  the  initialized  system  (sources  and  receivers)  is  reasonably  near  the  solution. 
This  is  the  case  when  the  source  positions  are  described  in  terms  of  x-y  locations  or 
bearings,  and  when  these  x-y  positions  are  accurately  initialized  using  GPS 
measurements. 

For  a  more  complicated  source  track  description,  such  as  the  4-parameter 
approach  discussed  in  Section  2.3,  the  error  surface  is  no  longer  quadratic  with  respect  to 
the  source  parameters,  and  so  the  Hessian-based  method  cannot  be  reliably  used.  In  such 
a  case,  a  steepest  descent  algorithm  is  employed,  the  update  to  the  unknown  parameters 
being  based  on  the  local  gradient  of  the  RMS  error.  This  algorithm  can  navigate  through 
a  more  complex  error  surface  in  search  of  the  optimal  solution. 
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Although  the  steepest  descent  algorithm  is  more  robust  against  complexities  in  the 
error  surface  and  against  poor  initialization,  the  inverse  Hessian  method  is  much  faster, 
and  so  is  used  where  possible.  The  inverse  Hessian  method  generally  requires 
10-25  iterations  to  converge,  while  the  steepest  descent  method  requires  at  least  50. 

2.3  ADS  AEL  ALGORITHM  VARIANTS 

This  section  describes  details  of  the  source  track  description,  the  relationship 
between  the  source/receiver  positions  and  the  time  delays,  and  the  means  of  initializing 
the  unknown  parameters.  The  algorithms  presented  here  can  be  organized  as  shown  in 
Fig.  2.3,  according  to  the  availability  of  accurate  initialization  and  the  distance  from  the 
receivers  to  the  source. 

The  availability  of  accurate  initialization  of  unknown  parameters  depends  on 
whether  dedicated  assets  that  can  provide  accurate  navigational  information  are  available 
for  AEL  runs,  or  whether  analysis  will  be  based  on  sources  of  opportunity,  such  as 
passing  merchant  shipping.  When  navigational  information  is  available,  source  positions 
can  be  accurately  initialized  directly  from  the  GPS  data.  In  these  cases,  GPS  data  are  also 
generally  available  for  the  deployment  platform,  allowing  for  an  accurate  estimate  of  the 
receiver  array's  absolute  position  and  line  of  bearing. 

When  accurate  navigational  information  is  not  available,  parameters  describing 
the  source  track  must  be  initialized  via  grid  search  methods,  obtaining  an  initialization 
that  minimizes  the  RMS  difference  between  the  measured  and  the  predicted  time  delays. 
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The  second  trait  distinguishing  AEL  algorithm  variants  from  one  another  is  the 
distance  from  the  receiver  array  to  the  acoustic  source  being  used  for  AEL.  One  option  is 
to  have  sources  in  the  near  field,  meaning  that  propagation  from  source  to  receiver  is  via 
direct  path.  This  lack  of  bottom  interactions  means  that  the  refracted  ray  paths  can  be 
accurately  modeled  even  when  no  information  regarding  bottom  structure  is  available. 

Generally,  detailed  bottom  information  is  not  available  in  the  early  phases  of  a 
deployment,  which  is  when  AEL  must  be  performed.  For  sources  in  direct  path  range, 
propagation  path  information  (including  travel  time)  can  be  obtained  using  only  the 
sound  speed  or  velocity  profile  (SVP)  in  the  water,  which  is  generally  available.  In  this 
case,  a  ray  trace  program  is  used  to  obtain  the  relationship  between  source-receiver 
separation  and  travel  time,  and  this  information  can  then  be  used  to  compute  modeled 
travel  time  differences  based  on  the  current  estimates  of  the  source  and  receiver  positions. 

AEL  analysis  can  also  be  based  on  far-field  sources.  There  is  no  clear  boundary 
between  far-field  and  mid-field  sources,  only  a  gradually  diminishing  error  resulting  from 
substituting  one  for  the  other.  Simply,  a  far-field  source  is  one  for  which  the  wavefront  of 
the  arriving  signal  can  be  treated  as  a  plane  wave.  For  a  source  that  is  15  km  distant,  the 
actual  wavefront  will  deviate  from  planar  by  about  2  m  across  a  500-m  aperture.  A 
deviation  of  1  m  would  require  a  contact  about  30  km  away.  The  far-field  AEL  algorithm 
would  represent  this  wavefront  as  planar,  and  so  induce  a  comparable  bow  in  the  receiver 
shape,  as  illustrated  in  Fig.  2.4.  Treating  a  wavefront  as  planar  greatly  simplifies  the 
algorithm,  resulting  in  improved  robustness  and  speed,  but  it  remains  important  to  estimate 
the  range  to  the  source  and,  if  necessary,  correct  the  overall  array  curvature. 
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As  shown  in  Fig.  2.3,  these  two  sets  of  options  allow  for  four  algorithm  choices, 
one  of  which  (far  field  with  navigation)  has  been  left  unimplemented  because  of 
irrelevance.  The  other  three  options  have  been  implemented  within  the  ADS  AEL  software 
suite. 


The  first  AEL  variant  is  a  near-field  source  with  known  navigation.  This  is  a 
variant  that  traditionally  has  shown  excellent  robustness  and  accuracy,  having  been  used 
successfully  on  dozens  of  data  sets,  both  by  ADS  and  other  programs.  This  algorithm 
requires  a  reasonably  accurate  estimate  of  the  sound  speed  profile  and  the  water  depth. 
Like  all  variants  discussed  here,  it  is  currently  implemented  to  treat  all  receiver  elements  as 
being  in  the  same  horizontal  plane,  although  this  variant  can  be  expanded  to  treat  the 
receiver  element  depths  as  unknowns. 

The  second  AEL  variant  is  a  far-field  source  with  unknown  navigation.  This 
algorithm  was  developed  specifically  for  the  ADS  program  in  response  to  the  need  for 
covert  AEL.  Each  source  position  is  fully  defined  by  its  bearing,  which  is  initialized  by  a 
simple  grid  search.  This  algorithm  has  proven  robust,  accurate,  and  fast,  with  the  primary 
source  of  error  being  the  incorrect  identification  of  a  mid-field  source  as  being  in  the  far 
field.  The  most  effective  way  to  rectify  incorrect  mid-field  identification  appears  to  be  by 
using  a  field  processor,  which  would  make  use  of  multi-node  holding  of  contacts  to  identify 
the  sources  of  opportunity  and  to  estimate  their  ranges. 

The  third  AEL  variant  is  a  near-field  source  with  unknown  navigation. 
Development  of  this  algorithm  was  motivated  not  only  by  the  need  for  covert  AEL  but 
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also  by  the  difficulty  of  identifying  true  far-field  contacts,  as  well  as  the  possible  lack  of 
high  signal-to-noise  ratio  (SNR)  far-field  contacts  in  highly  cluttered  environments. 
Near-field  sources  of  opportunity  are  available  in  abundance  in  most  environments; 
therefore  SNR  is  not  an  issue.  The  difficulty  of  this  method  is  that  too  little  information 
is  available  to  guarantee  algorithm  convergence  in  a  general  case.  To  limit  the  degrees  of 
freedom  of  the  system,  it  is  assumed  that  the  source  of  opportunity  has  passed  the  array 
on  a  constant-speed  straight-line  course.  Thus,  the  source  track  can  be  defined  by  four 
parameters,  in  this  case  being  the  time  of  closest  point  of  approach  (CPA),  the  bearing  to 
the  CPA  point,  the  CPA  range,  and  the  speed.  With  these  four  parameters  as  unknowns, 
along  with  the  x-y  location  of  each  receiver  element,  a  solution  can  usually  be  obtained. 
This  algorithm  has  not  proven  to  be  as  robust  and  reliable  as  the  other  two,  mainly  due  to 
the  difficulty  of  obtaining  an  accurate  enough  initialization  of  the  source  track,  and  so  it 
should  still  be  considered  to  be  in  a  developmental  stage. 

The  distinction  between  having  navigation  available  and  not  available  can 
generally  be  compared  to  the  difference  between  overt  and  covert  AEL.  Overt  AEL  is 
done  in  a  way  that  could  compromise  the  location  of  the  nodes,  such  as  driving  a  back- 
and-forth  pattern  around  each  node  with  a  surface  ship,  as  was  done  in  FET.  Covert  AEL 
indicates  the  intent  to  derive  array  configurations  without  doing  anything  that  would  give 
away  the  locations  of  the  nodes,  as  was  demonstrated  at  FET  using  the  far-field  no¬ 
navigation  algorithm.  There  are  exceptions  to  this,  such  as  using  radar  to  obtain 
navigational  information  on  near-field  sources  of  opportunity,  or  having  a  dedicated 
surface  asset  transit  through  the  field  in  a  fairly  linear  manner  but  along  a  path  that  brings 
it  past  some  or  all  of  the  nodes  in  the  field. 
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2.4  QUALITY  ASSESSMENT  OF  ALGORITHM  RESULTS 


As  the  AEL  process  proceeds,  various  diagnostic  displays  are  provided  to  allo'w 
for  an  assessment  of  the  validity  of  the  solution  obtained.  Some  of  these  are  listed  below. 

•  As  the  linetracker  runs,  it  displays  a  sequence  of  omnixomni 
correlagrams,  with  the  derived  time-T  relationship  superimposed.  This 
relationship  shows  if  the  linetracker  has  skipped  traces,  or  if  the  trace 
of  interest  has  faded  out. 

•  Scoring  of  the  time-x  points  is  shown,  which  allows  thresholding  to  be 
applied  to  remove  low  quality  data. 

•  As  AEL  iterates,  the  current  estimates  of  the  source  and  receiver 
positions  are  displayed.  For  the  near-field  navigation  method,  the 
distortion  of  the  source  track  should  be  minimal,  as  shown  in  Fig.  2.5. 
Significant  distortion  of  the  source  track,  as  for  example  in  Fig.  2.6, 
indicates  systematic  errors  such  as  depth  variation  of  the  receiver  array 
or  poor  initialization  of  sources  or  receivers.  Appendix  B  contains 
displays  for  all  events  of  the  source  track  upon  convergence. 

•  Upon  algorithm  convergence,  the  distances  between  adjacent  receiver 
elements  are  compared  to  the  planned  separations.  Since  the  length  of 
cable  between  elements  limits  their  true  separation,  significant  over- 
nominal  separation  is  generally  an  indication  of  problems  such  as 
incorrect  linetracking,  or  an  incorrect  environmental  description. 


14 


in 

o 

I 

o 

o 

CO 

< 


qwu  -  ui 


2 

+ 

c 

3 

ic 

zz 

o 

'TD 

C 

rt 

3 

o 

, 

C3 

a 

•  d 

s? 

'5) 

’C 

O 

00 

Itj 

r" 

.2 

5 

S 

00 

O 

c/: 

O 

Ou 

Cl. 

3 

> 

*5 

o 

o 

a:) 

cn 

0:^ 

8 

S 

> 

c 

o 


Sd 


(A 

s 

0 


o 

a 

0 

o  CJ 

wdS 

o 

.t  fl 

^  <u 

s  ^ 

S  ^ 

tN 

A  T-H 
V  ^ 
W  ^ 
^  0^ 
0  TJ 
O  O 

#  • 

V) 

<s 

U 


3d 

•  P4 


15 


Copyright®  1999 
The  University  of  Texas  at  Austin 
Applied  Research  Laboratories 
Reproduction  and  Redistribution  Prohibited 
Without  Prior  Express  Consent 


CD 

O 

I 

o 

o 

ob 

< 


C>  *-• 
l—  ZJ 

b  o 
d  t: 

T3  S 
C  ^ 

o  ^ 
^  a 
zz  c 

C3 

.E 

*5b*b 


CO 

C 

O 


00 

O 

D. 


(D 

U 


lx  © 

'O 


W 

C/3 


.  T3  ” 
©  ©  © 

^  i 

C5  C  ^ 


©  p«  ^ 

OX)  s  ^ 

^  m  pfi 
©  «  ^ 

Si  C3  OX) 
fi  fl 
O  i-»  *S 


o 

u  ^ 

,M5  C« 

1  S  ^ 

i-  «  «3 

O  ^  *s 

OX)  a  o 

«  P  cs 

(S  © 

5  T1  w) 

(M  OX) 

^  fi  5 
5«  .S  §: 

^  :S  ^ 

©  -73  IS 

S  g  2 

s  "  g 

^  t  ft 
fi  o  p 

»<  ©  © 

®  ©  « 

>..S  k( 

^  ^  O 

2  0^  ^ 
S  §  fa  S 

li-l^g 

.C  ft  o  c 

S  s  5  ^  §1 

5  «8  .2  to  « 
i  «  ^  C  £2 

g  ^  fTj 

s  5  2  2  *0 
©  «  2  ®  © 
««  nS  2 .2  z 

VO 

ri 

© 

u 

D 

OX) 


16 


Copyright(3)1999 
The  University  of  Texas  at  Austin 
Applied  Research  Laboratories 
Reproduction  and  Redistribution  Prohibited 
Without  Prior  Express  Consent 


•  Upon  algorithm  convergence,  the  RMS  error  between  the  modeled  and 
the  measured  time  delays  is  computed.  A  high  RMS  error  indicates 
either  premature  convergence  of  the  algorithm  or  noisy  but  otherwise 
correct  data. 

2.5  ARRAY  APERTURE  ESTIMATION  USING 
SOURCES  OF  OPPORTUNITY 

The  low  tension  deployment  of  the  FET  arrays  often  resulted  in  arrays  being 
relatively  straight  but  with  their  apertures  compressed.  When  these  arrays  were 
beamformed  with  a  nominal  (full  aperture)  array  description,  the  apparent  "endfire" 
beams  were  actually  oversteered,  resulting  in  reduced  sensitivity  of  those  beams.  This 
caused  white  space  to  appear  at  the  edges  of  the  correlagrams,  and  similarly  degraded  the 
narrowband  "endfire"  beams. 

Because  full  AEL  solutions  for  all  nodes  were  weeks  away  (due  mainly  to  RV7's 
busy  schedule),  a  means  was  needed  to  quickly  obtain  accurate  estimates  of  the  array 
aperture.  To  meet  this  need,  the  following  tool  was  developed  on-site.  A  series  of 
omnixomni  correlagrams  was  collected  (see  Fig.  2.7).  These  correlagrams  were  already 
being  computed  for  use  by  AEL.  Specific  omnidirectional  pairs  were  chosen  in  such  a 
way  that,  using  four  correlagrams,  the  arrival  time  at  hydrophone  1  could  be  related  to 
that  at  hydrophone  40  (of  the  40  LF  hydrophones)  for  use  in  estimating  the  entire 
aperture.  In  the  example  in  Fig.  2.7,  the  omnidirectional  pairs  are  1x13,  13x21,  21x31, 
and  31x40. 
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The  maximum  time  delay  (maximum  in  absolute  value)  of  acoustic  sources  was 
measured  on  each  correlagram.  It  was  assumed  that  some  source  of  opportunity  would 
cross  endfire  in  the  3  h  of  data  buffered  on  disk  at  any  one  time.  The  total  travel  time 
from  hydrophone  1  to  hydrophone  40  could  then  be  obtained  as  the  sum  of  the  four  travel 
times  extracted  from  these  omnixomni  correlagrams.  The  sound  speed  could  then  be 
used  to  convert  this  travel  time  to  an  aperture  estimate. 

This  process  could  then  be  repeated,  in  case  no  source  of  opportunity  had  actually 
reached  endfire  during  the  3-h  period.  In  fact,  because  the  entire  procedure  was  so  fast 
and  easy,  it  was  repeated  several  times  during  the  test,  generally  returning  very  consistent 
results. 


This  method  proved  to  be  very  accurate,  providing  an  even  better  estimate  of  the 
array  aperture  than  the  overt  AEL  method.  It  was  also  fast,  calculating  aperture  estimates 
for  all  16  nodes  in  a  few  hours. 
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3.  IMPLEMENTATION  IN  PAS 


The  AEL  algorithms  and  related  software  were  implemented  in  MATLAB  and  C, 
primarily  on  the  analysis  processor  (AP).  MATLAB  proved  a  very  useful  environment 
because  it  allowed  for  easy  code  modifications,  provided  convenient  graphics  tools,  and 
offered  easy  adaptation  to  new  analysis  tasks.  The  general  AP  configuration  and  data 
flow  are  summarized  in  Fig.  3.1. 

The  generic  multiprocessor  (GMP)  beamformer  produced  omnixomni 
correlagrams,  over  the  band  from  10  to  150  Hz,  for  a  set  of  omnidirectional  pairs  selected 
at  the  site  (to  account  for  dead  hydrophones,  etc.).  These  omnixomni  correlagrams  were 
made  available  in  shared  memory,  and  were  collected  by  a  disk  buffer  manager  resident 
on  the  AP,  and  stored  in  a  3-h  buffer  on  disk. 

A  MATLAB  program,  with  a  display  similar  to  that  in  Fig.  2.7,  imaged  the  data  in 
the  buffer.  When  data  useful  for  AEL  were  observed  on  a  particular  node,  the  buffered 
data  for  that  node  were  "grabbed"  into  a  permanent  disk  file  for  later  analysis. 

The  remaining  steps  in  the  AEL  process,  including  the  linetracker,  were 
implemented  in  MATLAB  and  C,  and  resided  exclusively  on  the  AP.  When  the  AEL 
algorithm  converged  to  a  solution,  and  that  solution  was  judged  acceptable  for  use  in  the 
remainder  of  PAS  (processing  automation  system),  a  configuration  file  was  generated  and 
copied  to  the  appropriate  beamformer. 
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4.  AEL  RESULTS  FOR  FET  AND  MNT 


The  AEL  results  for  all  FET  and  MNT  nodes  are  summarized  in  Figs.  4.1  to  4.5. 
In  these  figures,  the  dark  blue  and  cyan  (light  blue)  symbols  indicate  solutions  obtained 
with  the  overt  near-field  AEL  method,  using  RV7  or  RVl.  The  red  symbols  show  results 
from  the  far-field  AEL  algorithm,  which  used  distant  sources  of  opportunity.  Green 
symbols  show  results  from  the  element  calibration  (ELCAL)  analysis,  which  was  a 
completely  independent  array  shape  analysis  conducted  for  a  few  FET/MNT  arrays  by 
George  Shepard  at  BBN  Systems  and  Technologies.  Brown  symbols  designate  results 
from  the  near-field  covert  method,  which  is  still  under  development. 

For  all  methods,  the  following  convention  is  used  to  report  element  positions. 
The  jc  and  y  axes  are  in  the  horizontal  plane,  with  the  origin  of  coordinates  at  either 
hydrophone  1  or  at  the  lowest  numbered  normally  functioning  LF  hydrophone.  The  x- 
axis  passes  through  either  hydrophone  40  or  the  highest  numbered  normally  functioning 
LF  hydrophone.  The  x-y-z  axes  form  a  right-hand  coordinate  system,  with  the  z-axis 
pointing  into  the  ocean  floor.  So  that  this  document  may  be  unclassified,  the  x  and  y 
receiver  coordinates  are  normalized  by  the  constant  factor  fmin.  The  value  of  fmin  (or 
equivalently,  the  planned  array  aperture)  may  be  requested  from  PMS  407,  SPAWAR. 

All  AEL  results  make  use  of  the  aperture  estimates  collected  from  sources  of 
opportunity,  as  discussed  in  Section  2.5.  As  mentioned  there,  this  method  provided  the 
most  reliable  method  of  estimating  the  array  aperture,  and  all  results  obtained  from  other 
methods  were  scaled  accordingly. 
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In  general,  the  agreement  between  the  different  AEL  methods  is  quite  close.  In 
particular,  where  ELCAL  results  are  available,  they  generally  agreed  very  well  with  the 
corresponding  AEL  results.  There  was  sometimes  disagreement  between  the  far-field 
(covert)  and  the  near-field  (overt)  AEL  results,  generally  as  a  result  of  mistaking  a  mid- 
field  source  of  opportunity  for  one  in  the  far  field.  Performance  for  the  near-field  covert 
algorithm  was  less  satisfactory,  indicating  the  need  for  further  developmental  work. 

Figure  4.6  represents  tabulated  array  apertures  and  deviation  from  linear  for  all 
nodes,  along  with  a  note  indicating  what  array  configuration  was  used  in  the  PAS 
beamformer.  In  all  but  three  cases  the  array  was  sufficiently  straight,  taking  into 
account  such  factors  as  the  angle  between  the  two  subarrays,  to  use  a  linear  array 
description  in  PAS,  thereby  eliminating  the  problem  of  dealing  with  front  and  back 
beam  sets.  In  the  three  cases  where  the  arrays  were  not  straight  enough,  the  full  AEL 
solution  was  used  in  the  beamforming. 

The  lines  of  bearing  (LOBs)  collected  in  Fig.  4.6  note  the  position  of  hydrophone 
40  relative  to  hydrophone  1.  This  number  is  referred  to  as  the  direction  of  lay  (DOL)  in 
the  GMP  beamformer.  DOL  is  measured  in  units  of  degrees  east  of  north.  For  the  most 
part,  the  absolute  position  (latitude-longitude)  of  hydrophone  1  was  not  obtained  from 
AEL,  but  was  rather  an  input  to  the  AEL  process. 

In  Appendix  A,  the  best  estimate  of  x-y  hydrophone  positions  is  tabulated  for  each 
array,  using  the  convention  discussed  in  Section  2.5. 
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Node  %  full  aperture  Array  LOB  Node  Max  y-deviation  (m/finin)  Shape  used  in 

_ _ _ _ _ _  PAS  beamformer 

11  98.1  325®  MNTll  0.5  linear 

12  99.2  0®  MNT12 _ Og _ full  AEL 
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In  Appendix  B,  various  displays  are  shown  for  each  of  the  overt  near-field  AEL 
events.  These  displays  are  useful  for  determining  the  trustworthiness  of  the  results  for 
each  case.  Examples  of  the  omnixomni  correlagrams  are  shown,  which  demonstrate  the 
bearing  extent  achieved  by  the  source  (more  is  better)  and  the  SNR  (more  is  better).  It 
appeared  that  RVl  was  a  much  quieter  ship  than  RV7,  to  the  extent  that  RVl  was  nearly 
too  quiet  for  use  as  an  AEL  source.  RVl  was  used  for  the  AEL  runs  at  nodes  23,  41, 
and  42. 

Other  displays  in  Appendix  B  include  run  geometry,  both  from  the  navigation 
(top)  and  upon  algorithm  convergence  (bottom).  Results  suggest  that  good  agreement 
between  the  initial  and  the  converged  geometry  is  the  strongest  indicator  of  a  correct 
solution. 
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5. 


DISCUSSION  OF  AEL  RESULTS 


5.1  COMPARISON  OF  RESULTS  AND  ACCURACY  ESTIMATES 

Excluding  the  results  from  the  near-field  covert  method,  the  different  AEL 
algorithms  generally  provided  comparable  results.  One  exception  to  this  occurred  during 
use  of  the  far-field  covert  method,  when  a  midrange  source  of  opportunity  was  mistaken 
for  a  far-field  contact.  It  was  often  difficult  to  distinguish  between  these  two  cases  based 
solely  on  the  omnixomni  correlagram,  as  the  example  in  Fig.  5.1  illustrates.  This  source 
of  error  came  in  to  play  for  the  far-field  solutions  for  nodes  31,  32,  44,  and  perhaps  for 
MNT  1 1  as  well. 

As  discussed  earlier,  the  field  processor  can  provide  the  means  to  determine  the 
range  to  a  source  of  opportunity  being  considered  for  far-field  AEL.  If  the  source  is  not 
sufficiently  distant,  it  should  not  be  selected  for  use  in  AEL.  As  an  alternative  to 
discarding  such  a  source  of  opportunity,  the  data  from  that  source  could  be  processed  and 
then  a  correction  to  the  overall  array  curvature  could  be  applied  to  the  result,  based  on  the 
range  to  the  source  as  obtained  from  the  field  processor. 

For  cases  where  ELCAL  and  overt  AEL  results  were  both  available  (nodes  22  and 
MNT  12),  the  agreement  was  excellent.  Both  the  overall  shape  and  the  detailed 
variations  matched  well.  Since  agreement  between  such  qualitatively  different  methods 
cannot  be  coincidental,  these  two  cases  provide  the  best  opportunity  for  estimating  the 
accuracy  of  the  overt  AEL  results  in  general. 
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In  both  of  these  cases,  the  maximum  deviation  between  element  positions  found 
by  the  two  methods  is  about  0.5  m,  with  a  mean  deviation  of  about  0.3  m.  This  degree  of 
error  corresponds  to  10%  of  a  wavelength  at  300  to  500  Hz,  offering  precision  adequate 
for  the  ADS  application. 

These  two  nodes  (22  and  MNT  12)  represent  optimal  cases  for  the  AEL 
algorithm,  where  there  were  no  complicating  issues  (such  as  bottom  slopes),  and  where 
no  red  flags  were  present  in  the  operation  of  the  AEL  algorithm  (such  as  source  track 
distortion).  Some  arrays,  such  as  nodes  23,  31,  33,  and  44,  showed  indications  of  such 
difficulties,  and  so  should  be  assigned  a  somewhat  higher  uncertainty  (perhaps  1.0  m). 

5.2  AUTOMATION  OF  THE  AEL  PROCESS 

In  general,  during  FET  the  performance  reliability  and  accuracy  of  the  AEL  array 
shape  estimation  algorithm  were  excellent.  The  major  obstacle  remaining  before  it  is 
considered  "fleet  ready"  is  its  degree  of  automation.  That  is,  the  AEL  algorithm 
implemented  in  a  production  ADS  system  must  be  capable  of  returning  array  shape 
estimates  without  the  focused  attention  of  a  scientist.  Based  on  FET  experience,  the 
algorithm  is  definitely  not  at  that  point  yet. 

One  option  is  to  integrate  additional  automation  into  the  correlation-based  process 
that  was  used  at  FET.  The  process  was  hands-on  by  intent  during  FET  to  allow  a  careful 
evaluation  of  its  performance  and  accuracy.  Certainly,  additional  automation  could  be 
added  to  the  code.  For  instance,  a  more  automated  linetracking  algorithm  under 
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development  at  ARL:UT  has  shown  reasonable  reliability  when  tested  on  FET  data. 
Also,  algorithms  could  be  developed  to  recognize  the  occurrence  of  a  usable  CPA  event, 
based  on  bearing  rate  and  SNR.  Such  an  algorithm  is  in  the  developmental  stage  at 
ARL:UT. 

Already,  in  addition  to  operating  reliably  and  accurately,  the  entire  AEL 
algorithm  iterates  to  completion  without  requiring  analyst  interaction.  The  role  of  the 
analyst  is  to  examine  the  algorithm  output  and  assess  its  quality.  However,  the  primary 
quality  assessment  clues  (distortion  of  source  track  during  iteration,  over-nominal 
element  separations,  and  RMS  error  upon  convergence)  could  be  quantitatively  measured 
and  acceptability  thresholds  defined  for  each. 

In  this  scenario — using  these  acceptability  thresholds  as  screening  criteria — AEL 
would  be  running  out-of-sight  of  the  operator,  and  when  a  reasonable  solution  had  been 
obtained  for  a  particular  node,  an  alert  would  be  logged.  It  would  still  be  necessary  for  a 
somewhat  knowledgeable  human  operator  to  inspect  the  results,  which  in  this  case  would 
be  a  sample  omnixomni  correlagram  with  the  trace  selected,  the  quantitative  Q/A  clues 
with  thresholds  noted,  and  a  before-and-after  beam  correlagram.  Previous  solutions 
obtained  for  this  node  would  also  be  provided  for  comparison.  The  operator  would  then 
decide  whether  to  update  the  beamformer. 

In  addition  to  correlation-based  AEL,  alternative  algorithms  are  also  being 
considered.  Attention  has  been  focused  on  the  correlation-based  algorithm  described 
here  because  of  its  accuracy  and  reliability.  However,  other  array  shape  estimation 
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algorithms  have  also  shown  promise,  and  offer  better  prospects  for  automation.  For 
example,  a  signal-gain-based  algorithm  has  been  examined  at  ARL:UT,  and  has  shown 
good  performance  on  simulated  data.  Additional  research  on  this  and  other  alternative 
approaches  may  produce  an  algorithm  that  offers  acceptable  accuracy  with  a  higher 
degree  of  automation  than  correlation-based  AEL  can  provide. 

Finally,  the  array  aperture  estimation  method  that  uses  sources  of  opportunity 
could  probably  be  automated  fairly  easily  to  quickly  provide  good  aperture  estimates  with 
minimal  operator  interaction.  This  algorithm  can  be  viewed  as  the  first  step  in  the  overall 
AEL  process,  and  could  possibly  be  expanded  to  compute  all  of  the  receiver  x-positions 
(positions  parallel  to  the  array  axis).  Knowing  the  x-positions  would  make  it  easier  to 
obtain  the  y-positions  (perpendicular  to  array  axis)  using  other  methods. 

5.3  THREE-DIMENSIONAL  AEL 

It  became  apparent  during  FET  that  mission  planning  and  deployment 
considerations  prevent  placement  of  all  arrays  in  areas  where  the  bottom  is  flat.  This  is  a 
significant  factor  in  the  evaluation  of  the  AEL  algorithm,  since  for  a  few  arrays,  the 
known  depth  variation  over  the  aperture  of  the  array  accounted  for  some  discrepancies 
observed  in  the  AEL  results. 

This  finding  has  motivated  discussion  of  restructuring  AEL  to  solve  for  the  depth 
of  each  receiver,  rather  than  assuming  that  all  array  elements  are  at  the  same  depth. 
Earlier  versions  of  the  correlation-based  AEL  algorithm  have  solved  for  the  depth  of  each 
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receiver,  but  for  reasons  of  speed  and  robustness,  a  two-dimensional  approach  was 
implemented  in  PAS  during  FET.  Certainly,  it  is  possible  to  implement  a  full  three- 
dimensional  algorithm  for  ADS,  but  at  the  cost  of  some  complexity  and  computational 
load,  as  well  as  possible  reductions  in  the  robustness  of  the  solution. 

A  simpler  alternative  would  be  to  add  an  AEL  input  for  the  depth  profile  collected 
from  the  towed  deployment  vehicle  (TDV).  Comparison  of  the  depth  profile  obtained 
from  the  ELCAL  analysis  to  that  from  the  TDV  showed  good  agreement,  indicating  that 
the  TDV  depth  information  provided  a  reasonably  accurate  estimate.  With  this  as  an 
input  to  AEL,  a  solution  for  the  x  and  y  components  of  the  receiver  positions  would  be 
straightforward,  and  the  errors  associated  with  ignoring  the  depth  variations  would  be 
avoided. 
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Results  —  Best  Elstimates  of 
x-y  Hydrophone  Positions 
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BEST  ESTIMATES  OF 
x-y  HYDROPHONE  POSITIONS 


Element  positions  will  be  tabulated  for  each  array  deployed  during  FET.  These 
results  are  the  "best"  approximation,  generally  based  on  the  near-field  overt  AEL  analysis 
using  RVl  and  RV7,  although  the  ELCAL  results  are  used  for  node  MNT  11,  and  the 
near-field  covert  results  are  used  for  node  23. 

The  following  convention  is  used  to  report  element  positions.  The  x  and  y  axes 
are  in  the  horizontal  plane,  with  the  origin  of  coordinates  at  either  hydrophone  1  or  at  the 
lowest  numbered  non-dead  LF  hydrophone.  The  x-axis  passes  through  either  hydrophone 
40  or  the  highest  numbered  non-dead  LF  hydrophone.  The  x-y-z  axes  form  a  right-hand 
coordinate  system,  with  the  z-axis  pointing  into  the  ocean  floor. 
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NodeMNT  11 


Hydrophone  Number 


Planned  Positions 


Best  Solution  (ELCAL) 


1-40 


1-48 


1-52 


x(m/fmin) 


x(m/fniin) 


y(m/fmin) 


0.00 


x(m) 


y(m) 


1.67 


0.00 


0.00 


3.40 


-NOTLOCALIZED- 


11 


15 


5.28 


-NOTLOCALIZED- 


13 


17 


7.19 


-NOTLOCALIZED- 


14 


18 


9.22 


6.80 


-0.01 


15 


19 


11.36 


8.75 


-0.06 


16 


20 


13.60 


10.88 


-0.13 


17 


21 


15.96 


13.02 


-0.17 


10 


18 


22 


18.44 


15.35 


-0.19 


11 


19 


23 


21.04 


17.82 


-0.17 


12 


20 


24 


23.77 


-DEAD- 


13 


21 


25 


26.65 


22.89 


-0.22 


14 


22 


26 


29.67 


25.76 


-0.27 


15 

23 

27 

32.84 

28.73 

-0.25 

16 

24 

28 

36.17 

31.87 

-0.25 

17 

25 

29 

39.67 

35.15 

-0.21 

18 

26 

30 

43.35 

38.63 

-0.23 

19 

27 

31 

47.22 

42.25 

-0.19 

20 

28 

32 

51.29 

46.09 

-0.23 

21 

29 

33 

55.56 

50.16 

-0.16 

22 

30 

34 

57.22 

54.37 

-0.17 

23 

31 

35 

58.97 

56.03 

-0.13 

24 

32 

36 

60.81 

57.74 

-0.08 

25 

33 

37 

62.75 

59.57 

-0.11 

26 

34 

38 

64.78 

61.49 

-0.07 

27 

35 

39 

66.91 

63.47 

-0.05 

28 

36 

40 

69.16 

65.56 

-0.08 

29 

37 

41 

71.52 

67.76 

-0.08 

30 

38 

42 

73.99 

70.07 

-0.08 

31 

39 

43 

76.59 

72.42 

-0.03 

32 

40 

44 

79.33 

74.93 

-0.04 

33 

41 

45 

82.20 

77.57 

-0.06 

34 

42 

46 

85.22 

80.36 

-0.08 

35 

43 

47 

88.39 

83.29 

-0.12 

36 

44 

48 

91.73 

86.43 

-0.07 

37 

45 

49 

95.23 

-DEAD- 

38 

46 

50 

98.91 

93.16 

-0.03 

39 

47 

51 

102.78 

96.75 

-0.04 

40 

48 

52 

106.84 

100.49 

0.01 
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Node  MNT  12 


Hydrophone  Number  Planned  Positions 


x(ni/fmin) 

x(ni/fmin) 

0.00 

-DEAD- 

1.67 

0.00 

3.40 

1.59 

5.28 

3.27 

7.19 

5.01 

9.22 

6.90 

11.36 

8.88 

13.60 

11.07 

15.96 

13.26 

18.44 

15.63 

21.04 

18.08 

23.77 

20.72 

26.65 

23.39 

29.67 

26.26 

32.84 

29.29 

36.17 

32.46 

39.67 

-DEAD- 

43.35 

39.30 

47.22 

43.01 

51.29 

46.83 

55.56 

50.95 

57.22 

52.50 

58.97 

54.14 

60.81 

55.86 

62.75 

57.63 

64.78 

59.47 

66.91 

61.49 

69.16 

63.58 

71.52 

65.80 

73.99 

68.07 

76.59 

70.42 

79.33 

72.95 

82.20 

75.64 

85.22 

78,39 

88.39 

81.20 

91.73 

84.36 

95.23 

87.56 

98.91 

91.03 

102.78 

94.68 

106.84 

98.53 

Solution 
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Node  11 


Hydrophone  Number 

Planned  Positions 

1-40 

1-48 

1-52 

x(ni/fmin) 

1 

1 

1 

0.00 

2 

2 

2 

1.67 

3 

7 

7 

3.40 

4 

11 

15 

5.28 

5 

13 

17 

7.19 

6 

14 

18 

9.22 

7 

15 

19 

11.36 

8 

16 

20 

13.60 

9 

17 

21 

15.96 

10 

18 

22 

18.44 

11 

19 

23 

21.04 

12 

20 

24 

23.77 

13 

21 

25 

26.65 

14 

22 

26 

29.67 

15 

23 

27 

32.84 

16 

24 

28 

36.17 

17 

25 

29 

39.67 

18 

26 

30 

43.35 

19 

27 

31 

47.22 

20 

28 

32 

51.29 

21 

29 

33 

55.56 

22 

30 

34 

57.22 

23 

31 

35 

58.97 

24 

32 

36 

60.81 

25 

33 

37 

62.75 

26 

34 

38 

64.78 

27 

35 

39 

66.91 

28 

36 

40 

69.16 

29 

37 

41 

71.52 

30 

38 

42 

73.99 

31 

39 

43 

76.59 

32 

40 

44 

79.33 

33 

41 

45 

82.20 

34 

42 

46 

85.22 

35 

43 

47 

88.39 

36 

44 

48 

91.73 

37 

45 

49 

95.23 

38 

46 

50 

98.91 

39 

47 

51 

102.78 

40 

48 

52 

106.84 

Best  Solution 


x(ni/fmin) 

■R9SI 

0.00 

0.00 

1.61 

-0.12 

3.25 

-0.12 

-DEAD- 

6.81 

0.04 

-DEAD- 

-DEAD- 

13.11 

0.14 

15.44 

0.12 

17.88 

0.13 

20.37 

0.08 

23.09 

0.00 

25.94 

-0.01 

28.86 

0.04 

32.03 

0.09 

35.32 

0.10 

38.79 

42.44 

46.27 

0.06 

50.29 

0.03 

54.54 

0.07 

56.15 

0.09 

57.91 

0.07 

59.75 

0.05 

61.67 

0.04 

63.64 

0.05 

65.74 

0.08 

67.91 

0.10 

70.24 

0.12 

72.69 

0.09 

75.21 

0.09 

77.88 

0.11 

80.76 

0.09 

83.74 

0.08 

86.85 

0.10 

90.20 

0.10 

93.73 

0.06 

97.33 

0.00 

101.14 

-0.02 

104.77 

0.00 

45 


Hydrophone  Number 

Planned  Positions 

1-40 

1-48 

1-52 

x(ni/fmin) 

1 

1 

1 

0.00 

2 

2 

2 

1.67 

3 

7 

7 

3.40 

4 

11 

15 

5.28 

5 

13 

17 

7.19 

6 

14 

18 

9.22 

7 

15 

19 

11.36 

8 

16 

20 

13.60 

9 

17 

21 

15.96 

10 

18 

22 

18.44 

11 

19 

23 

21.04 

12 

20 

24 

23.77 

13 

21 

25 

26.65 

14 

22 

26 

29.67 

15 

23 

27 

32.84 

16 

24 

28 

36.17 

17 

25 

29 

39.67 

18 

26 

30 

43.35 

19 

27 

31 

47.22 

20 

28 

32 

51.29 

21 

29 

33 

55.56 

22 

30 

34 

57.22 

23 

31 

35 

58.97 

24 

32 

36 

60.81 

25 

33 

37 

62.75 

26 

34 

38 

64.78 

27 

35 

39 

66.91 

28 

36 

40 

69.16 

29 

37 

41 

71.52 

30 

38 

42 

73.99 

31 

39 

43 

76.59 

32 

40 

44 

79.33 

33 

41 

45 

82.20 

34 

42 

46 

85.22 

35 

43 

47 

88.39 

36 

44 

48 

91.73 

37 

45 

49 

95.23 

38 

46 

50 

98.91 

39 

47 

51 

102.78 

Best  Solution 


x(ni/fmin) 

■msssfi 

0.00 

0.00 

1.60 

-0.06 

3.18 

-0.10 

4.82 

-0.02 

-DEAD-- 

-DEAD- 

10.58 

0.25 

12.74 

0.29 

15.11 

0.33 

17.61 

0.31 

-DEAD- 

22.66 

0.30 

25.47 

0.29 

28.51 

0.22 

31.61 

0.19 

34.87 

0.21 

38.37 

0.17 

42.06 

0.17 

45.92 

0.15 

50.00 

0.13 

54.31 

0.12 

55.91 

0.07 

57.71 

0.00 

59.54 

-0.02 

61.46 

-0.01 

63.51 

0.00 

65.63 

0.02 

67.85 

-0.02 

-DEAD- 

72.60 

-0.05 

75.29 

0.00 

77.96 

-0.01 

80.91 

-0.03 

83.95 

-0.09 

87.27 

-0.09 

90.58 

-0.07 

94.10 

-0.08 

97.91 

-0.03 

-DEAD- 

Node  13 


Hydro] 

phone  > 

umber 

Planned  Positions 

Best  Solution 

1-40 

1-48 

1-52 

x(ni/fnun) 

x(in/fmin) 

y(ni/fmin) 

1 

1 

1 

0.00 

0.00 

0.00 

2 

2 

2 

1.67 

1.61 

0.10 

3 

7 

7 

3.40 

3.12 

0.01 

4 

11 

15 

5.28 

4.72 

-0.15 

5 

13 

17 

7.19 

6.31 

0.25 

6 

14 

18 

9.22 

8.12 

0.21 

7 

15 

19 

11.36 

9.97 

0.30 

8 

16 

20 

13.60 

-DEAD- 

9 

17 

21 

15.96 

14.13 

0.51 

10 

18 

22 

18.44 

16.39 

0.38 

11 

19 

23 

21.04 

-DEAD- 

12 

20 

24 

23.77 

21.41 

0.35 

13 

21 

25 

26.65 

24.00 

0.07 

14 

22 

26 

29.67 

-DEAD- 

15 

23 

27 

32.84 

29.91 

0.02 

16 

24 

28 

36.17 

33.09 

0.09 

17 

25 

29 

39.67 

36.51 

0.08 

18 

26 

30 

43.35 

40.03 

0.19 

19 

27 

31 

47.22 

43.74 

0.28 

20 

28 

32 

51.29 

47.60 

0.02 

21 

29 

33 

55.56 

51.80 

0.21 

22 

30 

34 

57.22 

53.42 

0.10 

23 

31 

35 

58.97 

55.15 

0.17 

24 

32 

36 

60.81 

56.98 

0.24 

25 

33 

37 

62.75 

59.04 

0.55 

26 

34 

38 

64.78 

-DEAD- 

27 

35 

39 

66.91 

-DEAD- 

28 

36 

40 

69.16 

64.65 

-0.11 

29 

37 

41 

71.52 

66.82 

0.10 

30 

38 

42 

73.99 

69.03 

0.08 

31 

39 

43 

76.59 

71.29 

-0.05 

32 

40 

44 

79.33 

-NOT  LOCALIZED- 

33 

41 

45 

82.20 

-NOT  LOCALIZED- 

34 

42 

46 

85.22 

78.80 

0.05 

35 

43 

47 

88.39 

81.43 

0.13 

36 

44 

48 

91.73 

84.46 

0.18 

45 

49 

95.23 

87.72 

0.02 

38 

46 

50 

98.91 

91.04 

-0.04 

39 

47 

51 

102.78 

94.49 

0.16 

48 

52 

106.84 

98.15 

0.00 

•7 


Node  14 


Hydrophone  Number 

Planned  Positions 

Best  Solution 

1-40 

1-48 

1-52 

x(ni/fmin) 

x(m/fmin) 

1 

1 

1 

0.00 

0.00 

0.00 

2 

2 

2 

1.67 

1.56 

-0.21 

3 

7 

7 

3.40 

3.09 

-0.23 

4 

11 

15 

5.28 

4.61 

-0.13 

5 

13 

17 

7.19 

6.19 

0.12 

6 

14 

18 

9.22 

7.99 

0.46 

7 

15 

19 

11.36 

9.86 

0.39 

8 

16 

20 

13.60 

11.92 

0.49 

9 

17 

21 

15.96 

14.08 

0.30 

10 

18 

22 

18.44 

16.35 

0.27 

11 

19 

23 

21.04 

18.68 

0.33 

12 

20 

24 

23.77 

21.26 

0.49 

13 

21 

25 

26.65 

--DEAD- 

14 

22 

26 

29.67 

26.50 

15 

23 

27 

32.84 

29.37 

1.07 

16 

24 

28 

36.17 

32.38 

1.22 

17 

25 

29 

39.67 

35.68 

1.26 

18 

26 

30 

43.35 

39.04 

1.34 

19 

27 

31 

47.22 

42.62 

1.32 

20 

28 

32 

51.29 

46.36 

1.34 

21 

29 

33 

55.56 

50.15 

1.21 

22 

30 

34 

57.22 

51.52 

1.11 

23 

31 

35 

58.97 

53.11 

1.14 

24 

32 

36 

60.81 

54.66 

1.11 

25 

33 

37 

62.75 

56.24 

1.15 

26 

34 

38 

64.78 

57.97 

1.05 

27 

35 

39 

66.91 

59.68 

1.15 

28 

36 

40 

69.16 

61.53 

1.19 

29 

37 

41 

71.52 

63.46 

1.37 

30 

38 

42 

73.99 

65.52 

1.49 

31 

39 

43 

76.59 

-DEAD- 

32 

40 

44 

79.33 

-DEAD- 

33 

41 

45 

82.20 

72.35 

1.08 

34 

42 

46 

85.22 

74.96 

0.97 

35 

43 

47 

88.39 

77.61 

0.81 

36 

44 

48 

91.73 

80,47 

0.64 

37 

45 

49 

95.23 

83.47 

0.38 

38 

46 

50 

98.91 

86.58 

0.21 

39 

47 

51 

102.78 

89.89 

0.13 

40 

48 

52 

106.84 

93.34 

0.00 

48 


Node  21 


Hydrophone 


1-40  1-48 


Number  Planned  Positions 


1-52 

x(ni/fmin) 

x(ni/fmin) 

1 

0.00 

-DEAD- 

2 

1.67 

0.00 

7 

3.40 

1.56 

15 

5.28 

3.39 

17 

7.19 

5.22 

18 

9.22 

7.12 

19 

11.36 

-DEAD- 

20 

13.60 

11.52 

21 

15.96 

13.83 

22 

18.44 

16.47 

23 

21.04 

-DEAD- 

24 

23.77 

21.72 

25 

26.65 

24.50 

26 

29.67 

27.58 

27 

32.84 

30.70 

28 

36.17 

34.10 

29 

39.67 

37.56 

30 

43.35 

41.34 

31 

47.22 

45.31 

32 

51.29 

49.38 

33 

55.56 

53.71 

34 

57.22 

55.38 

35 

58.97 

57.20 

36 

60.81 

59.17 

37 

62.75 

61.03 

38 

64.78 

62.84 

39 

66.91 

65.08 

40 

69.16 

67.17 

41 

71.52 

69.60 

42 

73.99 

72.06 

43 

76.59 

74.70 

44 

79.33 

77.45 

45 

82.20 

80.37 

46 

85.22 

83.39 

47 

88.39 

86.63 

48 

91.73 

89.98 

49 

95.23 

93.52 

50 

98.91 

97.19 

51 

102.78 

101.09 

52 

106.84 

105.17 

Best  Solution 


'(m/fmin) 


49 


Node  22 


Hydrophone  Number 

Planned  Positions 

Best  Solution 

1-40 

1-48 

1-52 

x(m/fmin) 

x(ni/fmin) 

IBBSSI 

1 

1 

1 

0.00 

0.00 

0.00 

2 

2 

2 

1.67 

-NOT  LOCALIZED- 

3 

7 

7 

3.40 

-NOT  LOCALIZED- 

4 

11 

15 

5.28 

-NOT  LOCALIZE-D- 

5 

13 

17 

7.19 

6.67 

0,00 

6 

14 

18 

9.22 

8.63 

0.00 

7 

15 

19 

11.36 

10.68 

-0.03 

8 

16 

20 

13.60 

12.83 

-0.10 

9 

17 

21 

15.96 

15.13 

-0.11 

10 

18 

22 

18.44 

17.52 

-0.18 

11 

19 

23 

21.04 

20.02 

-0.18 

12 

20 

24 

23.77 

22.72 

-0.21 

13 

21 

25 

26.65 

25.56 

-0.24 

14 

22 

26 

29.67 

28.51 

-0.29 

15 

23 

27 

32.84 

31.61 

-0.31 

16 

24 

28 

36.17 

34.88 

-0.41 

17 

25 

29 

39.67 

38.38 

-0.43 

18 

26 

30 

43.35 

42.03 

-0.49 

19 

27 

31 

47.22 

45.86 

-0.48 

20 

28 

32 

51.29 

49.94 

-0.48 

21 

29 

33 

55.56 

54,21 

-0.49 

22 

30 

34 

57.22 

55.83 

-0.48 

23 

31 

35 

58.97 

57.63 

-0.45 

24 

32 

36 

60.81 

59.45 

-0.42 

25 

33 

37 

62.75 

61.40 

-0.37 

26 

34 

38 

64.78 

63.44 

-0.34 

27 

35 

39 

66.91 

65.48 

-0.32 

28 

36 

40 

69.16 

67.76 

-0.36 

29 

37 

41 

71.52 

70.07 

-0.40 

30 

38 

42 

73.99 

72.48 

-0.45 

31 

39 

43 

76.59 

-DEAD- 

32 

40 

44 

79.33 

77.79 

-0.51 

33 

41 

45 

82.20 

80.62 

-0.53 

34 

42 

46 

85.22 

83.62 

-0.50 

35 

43 

47 

88.39 

86.88 

-0.48 

36 

44 

48 

91.73 

90.22 

-0.40 

37 

45 

49 

95.23 

93.72 

-0.35 

38 

46 

50 

98.91 

97.45 

-0.24 

39 

47 

51 

102.78 

101.32 

-0.09 

40 

48 

52 

106.84 

105.41 

0.00 
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Node  23 


Hydrophone  Number 

Planned  Posit! 

1-40 

1-48 

1-52 

x(m/fmin) 

1 

1 

1 

2 

2 

2 

1.67 

3 

7 

7 

3.40 

4 

11 

15 

5.28 

5 

13 

17 

7.19 

6 

18 

9.22 

7 

19 

11.36 

8 

16 

20 

9 

17 

21 

15.96 

10 

18 

22 

18.44 

11 

19 

23 

21.04 

12 

24 

23.77 

13 

21 

25 

26.65 

14 

22 

26 

29.67 

15 

23 

27 

32.84 

16 

24 

28 

36.17 

17 

25 

29 

39.67 

18 

26 

30 

43.35 

19 

27 

31 

47.22 

20 

28 

32 

51.29 

21 

29 

33 

55.56 

22 

30 

34 

57.22 

23 

35 

58.97 

24 

32 

36 

60.81 

25 

33 

37 

62.75 

26 

34 

38 

64.78 

27 

35 

39 

66.91 

28 

36 

40 

69.16 

29 

37 

41 

71.52 

30 

38 

42 

73.99 

31 

43 

76.59 

32 

44 

79.33 

33 

41 

45 

82.20 

34 

42 

46 

85.22 

35 

43 

47 

88.39 

36 

44 

48 

91.73 

37 

45 

49 

95.23 

38 

46 

98.91 

39 

47 

51 

40 

48 

52 

llllllllllllilllll^^ 

51 


Node  24 


52 


Node  32 


54 


Node  41 


Planned  Positions 


:  Solution 


1-40 

1-48 

1-52 

x(ni/fmin) 

x(m/fmin) 

y(ni/fniin) 

1 

1 

1 

0.00 

0.00 

0.00 

2 

2 

2 

1.67 

1.95 

-0.26 

3 

7 

7 

3.40 

3.25 

0.00 

7.19 


9.22 


11.36 


13.60 


15.96 


18.44 


21.04 


23.77 


26.65 


29.67 


32.84 


36.17 


39.67 


43.35 


47.22 


51.29 


55.56 


57.22 


58.97 


60.81 


62.75 


64.78 


66.91 


69.16 


71.52 


73.99 


76.59 


79.33 


82.20 


85.22 


88.39 


91.73 


95.23 


98.91 


102.78 


106.84 


6.88 


8.87 


11.03 


13.28 


15.67 


17.85 


20.42 


23.28 


26.17 


28.91 


32.21 


35.66 


39.16 


42.99 


46.91 


50.99 


55.15 


57.01 


58.81 


60.72 


62.79 


64.45 


66.59 


68.93 


71.20 


73.71 


76.35 


78.80 


81.67 


84.65 


87.76 


91.39 


95.00 


98.67 


102.71 


106.80 


55 


56 


Node  43 


Hydrophone  Number 

Planned  Positions 

Best  Solution 

1-40 

1-48 

1-52 

x(m/fmin) 

x(m/fmin) 

y(m/fmin) 

1 

1 

1 

0.00 

0.00 

0.00 

2 

2 

2 

1.67 

1.59 

0.12 

3 

7 

7 

3.40 

3.17 

0.37 

4 

11 

15 

5.28 

4.83 

0.42 

5 

13 

17 

7.19 

6.62 

0.29 

6 

14 

18 

9.22 

8.47 

0.29 

7 

15 

19 

11.36 

10.53 

8 

16 

20 

13.60 

0.38 

9 

17 

21 

15.96 

14.95 

0.43 

10 

18 

22 

18.44 

17.45 

0.45 

11 

19 

23 

21.04 

20.19 

0.49 

12 

20 

24 

23.77 

22.86 

13 

21 

25 

26.65 

25.52 

0.47 

14 

22 

26 

29.67 

28.49 

0.50 

15 

23 

imiiii 

32.84 

31.74 

0.54 

16 

24 

28 

36.17 

35.09 

0.57 

17 

25 

29 

39.67 

38.44 

nmg^i[iiiim 

18 

26 

30 

43.35 

42.10 

0.63 

19 

27 

31 

47.22 

45.85 

0.68 

20 

28 

32 

51.29 

49.89 

0.71 

21 

29 

33 

55.56 

54.09 

0.79 

22 

30 

34 

57.22 

55.75 

23 

31 

35 

58.97 

57.42 

0.82 

24 

32 

36 

60.81 

59.18 

0.81 

25 

33 

37 

62.75 

61.07 

0.87 

26 

34 

38 

64.78 

63.24 

0.85 

27 

35 

39 

66.91 

65.30 

0.89 

28 

36 

40 

69.16 

67.44 

0.87 

29 

37 

41 

71.52 

69.62 

0.87 

30 

38 

42 

73.99 

71.94 

0.90 

31 

39 

43 

76.59 

74.38 

0.87 

32 

40 

44 

79.33 

77.16 

0.77 

33 

41 

45 

82.20 

mewm 

34 

42 

46 

85.22 

82.59 

0.71 

35 

43 

47 

88.39 

85.46 

0.65 

36 

44 

48 

91.73 

88.61 

0.55 

37 

45 

49 

95.23 

92.02 

0.43 

46 

50 

98.91 

95.50 

0.35 

39 

47 

51 

102.78 

99.27 

0.21 

40 

48 

52 

106.84 

0.00 
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Node  44 


Hydrophone  Number 

Planned  Positions 

1-40 

1-48 

1-52 

x(m/fmin) 

1 

1 

1 

0.00 

2 

2 

2 

1.67 

3 

7 

7 

3.40 

4 

11 

15 

5.28 

5 

13 

17 

7.19 

6 

14 

18 

9.22 

7 

15 

19 

11.36 

8 

16 

20 

13.60 

9 

17 

21 

15.96 

10 

18 

22 

18.44 

11 

19 

23 

21.04 

12 

20 

24 

23.77 

13 

21 

25 

26.65 

14 

22 

26 

29.67 

15 

23 

27 

32.84 

16 

24 

28 

36.17 

17 

25 

29 

39.67 

18 

26 

30 

43.35 

19 

27 

31 

47.22 

20 

28 

32 

51.29 

21 

29 

33 

55.56 

22 

30 

34 

57.22 

23 

31 

35 

58.97 

24 

32 

36 

60.81 

25 

33 

37 

62.75 

26 

34 

38 

64.78 

27 

35 

39 

66.91 

28 

36 

40 

69.16 

29 

37 

41 

71.52 

1 

30 

38 

42 

73.99 

31 

39 

43 

76.59 

32 

40 

44 

79.33 

33 

41 

45 

82.20 

34 

42 

46 

85.22 

35 

43 

47 

88.39 

36 

44 

48 

91.73 

37 

45 

49 

95.23 

38 

46 

50 

98.91 

39 

47 

51 

102.78 

40 

48 

52 

106.84 

Best  Solution 


x(m/fmin) 

y(m/fmin) 

0.00 

0.00 

1.64 

0.17 

-DEAD- 

-DEAD- 

6.92 

0.29 

8.91 

0.23 

11.00 

0.17 

13.18 

0.10 

15.46 

0.09 

17.72 

0.06 

20.19 

-0.02 

22.93 

-0.06 

25.73 

-0.16 

28.68 

-0.17 

31.86 

-0.14 

35.16 

0.00 

38.63 

0.10 

42.36 

0.22 

46.14 

0.32 

50.11 

0.37 

54.25 

0.38 

55.82 

0.42 

57.55 

0.43 

59.40 

0.45 

61.35 

0.44 

63.27 

0.44 

65.33 

0.43 

67.49 

0.41 

0.43 

72.24 

0.42 

74.67 

0.39 

-DEAD- 

80.15 

0.27 

83.13 

0.20 

86.19 

0.15 

89.59 

0.10 

93.10 

0.14 

96.74 

0.12 

100.54 

0.09 

104.60 

0.00 
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APPENDIX  B 

Displays  of  Source  Tracks  upon  Convergence 
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h/min 


20/00 


Figure  B.l:  Nodell.  Omnixomni  correlagram  for  near-field 

contact  during  AEL  event,  with  processed  track  highlighted. 
Hydrophones  1  and  13,  day/h/min:  084/20/00  to  084/21/40 


AS-00-16 
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^  1 


Longitude  -  deg 

(a)  Navigation  geometry  during  AEL  event,  where  a  green  "+"  indicates  time  processed,  a  black  "+" 
indicates  time  not  processed,  and  a  red  "o"  indicates  receiver  location. 


m  -  east 

(b)  Source  and  receiver  positions  upon  convergence,  where  green  indicates  source  location,  red  indicates 
receiver  location,  +  indicates  estimated  location,  and  "o"  indicates  original  location. 

Figure  B.2:  Node  11.  Navigation  geometry  during  AEL  event, 

and  source  track  upon  convergence.  ao  nn  17 

Day/h/min:  084/20/00 

Copyright©  1999 
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h/min 


Time  delay  -  ms 


Figure  B.3:  Node  12.  Omnixomni  correlagram  for  near-field  contact  during 
AEL  event,  with  processed  track  highlighted. 

Hydrophones  1  and  13,  day/h/min:  085/05/30  to  085/08/00 
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33.38 


Longitude  -  deg 

(a)  Navigation  geometry  during  AEL  event,  where  a  green  ”+"  indicates  time  processed,  a  black  "+” 
indicates  time  not  processed,  and  a  red  "o"  indicates  receiver  location. 


(b)  Source  and  receiver  positions  upon  convergence,  where  green  indicates  source  location,  red  indicates 
receiver  location,  "+"  indicates  estimated  location,  and  "o"  indicates  original  location. 

Figure  B.4;  Node  12.  Geometry  during  AEL  event  (top)  and  source  track 

upon  convergence  (bottom).  aq  nn  i  q 

Day/h/min:  085/05/30  Ab-ou- 1 9 
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Figure  B.5:  Node  21.  Omnixomni  correlagram  for  near-Held  contact  during  AEL 
event,  with  processed  track  highlighted. 

Hydrophones  13  and  20,  day/h/min:  090/14/40  to  090/16/30 
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(a)  Navigation  geometry  during  AEL  event,  where  a  green  ”+"  indicates  time  processed,  a  black  "+" 
indicates  time  not  processed,  and  a  red  "o"  indicates  receiver  location. 
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(b)  Source  and  receiver  positions  upon  convergence,  where  green  indicates  source  location,  red  indicates 
receiver  location, "+"  indicates  estimated  location,  and  "o"  indicates  original  location. 


Figure  B,6:  Node  21.  Geometry  during  AEL  event  (top)  and  source 
track  upon  convergence  (bottom), 
day/h/min:  090/14/40 
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Figure  B.7:  Node  22.  Omnixomni  correlagram  for  near-field  contact  during  AEL 
event,  with  processed  track  highlighted. 

Hydrophones  2  and  13,  day/h/min:  107/07/00  to  107/10/00 


AS-00-22 
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(a)  Navigation  geometry  during  AEL  event,  \vhere  a  green  indicates  time  processed,  a  black 
indicates  time  not  processed,  and  a  red  "o”  indicates  receiver  location. 


(b)  Source  and  receiver  positions  upon  convergence,  where  green  indicates  source  location,  red  indicates 
receiver  location, indicates  estimated  location,  and  ”0"  indicates  original  location. 

Figure  B.8:  Node  22.  Geometry  during  AEL  event  (top)  and  source 

track  upon  convergence  (bottom).  AS- 00-23 

Day/h/min:  107/07/00 
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(b)  Source  receiver  positions  upon  convergence,  where  green  indicates  source  location,  red  indicates 
receiver  location,  +  indicates  estimated  location,  and  "o”  indicates  original  location. 


Figure  B.IO:  Node  23.  Geometry  during  AEL  event  (top)  and  source 
track  upon  convergence  (bottom). 

Day/h/min:  130/08/00 
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Figure  B.ll:  Node  24.  Omnixomni  correlagram  for  near-field  contact  during  AEL 
event,  with  processed  track  highlighted. 

Hydrophones  1  and  13,  day/h/min:  090/07/00  to  090/10/00 


AS-00-26 
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(a)  Navigation  geometry  during  AEL  event,  where  a  green  "+"  indicates  time  processed,  a  black "+" 
indicates  time  not  processed,  and  a  red  "o"  indicates  receiver  location. 
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(b)  Source  and  receiver  positions  upon  convergence,  where  green  indicates  source  location,  red  indicates 
receiver  location, "+"  indicates  estimated  location,  and  "o"  indicates  original  location. 


Figure  B.12:  Node  24.  Geometry  during  AEL  event  (top)  and  source 
track  upon  convergence  (bottom). 

Day/h/min:  190/07/00 
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Figure  B.13:  Node  31.  Omnixomni  correlagram  for  near-Held  contact  during  AEL 
event,  with  processed  track  highlighted. 

Hydrophones  1  and  13,  day/h/min:  087/05/30  to  087/08/00 


AS-00-28 
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(a)  Navigation  geometry  during  AEL  event,  where  a  green  "+”  indicates  time  processed,  a  black "+" 
indicates  time  not  processed,  and  a  red  "o"  indicates  receiver  location. 


(b)  Source  and  receiver  positions  upon  convergence,  where  green  indicates  source  location,  red  indicates 
receiver  location,  +  indicates  estimated  location,  and  "o"  indicates  original  location. 


Figure  B.14:  Node  31.  Geometry  during  AEL  event  (top)  and  source 
track  upon  convergence  (bottom). 

Day/h/min:  087/05/30 
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Figure  B.15:  Node  32.  Omnixomni  correlagram  for  near-field  contact  during  AEL 
event,  with  processed  track  highlighted. 

Hydrophones  1  and  13,  day/h/min:  087/07/00  to  087/09/30 


AS-00-30 
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(a)  Navigation  geometry  during  AEL  event,  where  a  green  "+”  indicates  time  processed,  and  a  red  "o" 
indicates  receiver  location. 
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(b)  Source  and  receiver  positions  upon  convergence,  where  green  indicates  source  location,  red  indicates 
receiver  location, "+"  indicates  estimated  location,  and  "o"  indicates  original  location. 


Figure  B.16:  Node  32.  Geometry  during  AEL  event  (top)  and  source 
track  upon  convergence  (bottom). 

Day/h/min:  087/07/00 
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Figure  B.17:  Node  33.  Omnixomni  correlagram  for  near-fleld  contact  during  AEL 
event,  with  processed  track  highlighted. 

Hydrophones  1  and  13,  day/h/min:  087/09/00  to  087/11/30 


AS-00-32 
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Figure  B.18:  Node  33.  Geometry  during  AEL  event  (top)  and  source 
track  upon  convergence  (bottom). 

Day/h/min:  087/09/00  AS-00-33 
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Figure  B.19:  Node  34.  Omni  correlagram  during  AEL  event,  with 

processed  track  highlighted.  Only  one  pass  is  illustrated,  but  three 
passes  were  processed  to  provide  the  information  for  Fig.  B.20. 
Hydrophones  1  and  13,  day/h/min:  094/19/30  to  094/20/10 


AS-00-34 
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(a)  Navigation  geometry  during  AEL  event,  where  a  green  "+"  indicates  time  processed,  a  black  "+" 
indicates  time  not  processed,  and  a  red  "o"  indicates  receiver  location. 


(b)  Source  ^d  receiver  positions  upon  convergence,  where  green  indicates  source  location,  red  indicates 
receiver  location,  +  indicates  estimated  location,  and  "o"  indicates  original  location. 


Figure  B.20:  Node  34.  Geometry  during  AEL  event  (top)  and  source 

track  upon  convergence  (bottom).  AS-00-35 

Day/h/min:  094/18/20 
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Figure  B.21;  Node  41.  Omnixomni  correlagram  for  near-fleld  contact  during  AEL 
event,  with  processed  track  highlighted. 

Hydrophones  1  and  13,  day/h/min:  130/10/00  to  130/11/10 


AS-00-36 
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(a)  Navigation  geometry  during  AEL  event,  where  a  green  "+"  indicates  time  processed,  a  black 
indicates  time  not  processed,  and  a  red  "o"  indicates  receiver  location. 


(b)  Source  and  receiver  positions  upon  convergence,  where  green  indicates  source  location,  red  indicates 
receiver  location, indicates  estimated  location,  and  "o"  indicates  original  location. 


Figure  B.22:  Node  41.  Geometry  during  AEL  event  (top)  and  source 
track  upon  convergence  (bottom). 

Day/h/min:  130/10/00 
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Figure  B.23;  Node  42.  Omnixomni  correlagram  for  near-fleld  contact  diu-ing  AEL 
event,  with  processed  track  highlighted. 

Hydrophones  1  and  13,  day/h/min:  130/13/00  to  130/14/50 


AS-00-38 
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Figure  B.24:  Node  42.  Navigation  geometry  during  AEL  event,  where  a  green  ’*+" 
indicates  time  processed,  a  black  indicates  time  not  processed, 
and  a  red  "o"  indicates  receiver  location. 

Day/h/min:  130/13/00 
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Figure  B.25:  Node  43.  Omnixomni  correlagram  for  near-field  contact  during  AEL 
event,  with  processed  track  highlighted. 

Hydrophones  1  and  13,  day/h/min:  101/20/40  to  101/23/30 


AS-00-40 
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(a)  Navigation  geometry  during  AEL  event,  where  a  green  indicates  time  processed,  a  black  ”+' 
indicates  time  not  processed,  and  a  red  "o”  indicates  receiver  location. 


(b)  Source  and  receiver  positions  upon  convergence,  where  green  indicates  source  location,  red  indicates 
receiver  location,  +  indicates  estimated  location,  and  ”0”  indicates  original  location. 

Figure  B.26:  Node  43.  Geometry  during  AEL  event  (top)  and  source 

track  upon  convergence  (bottom).  Ac;-nn-4i 

Day/h/min:  101/20/40 
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Figure  B.27:  Node  44.  Omnixomni  correlagram  for  near-field  contact  during  aEL 
event,  with  processed  track  highlighted. 

Hydrophones  1  and  13,  day/h/min:  084/15/40  to  084/17/30 


AS-00-42 
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convergence,  where  green  indicates  source  location,  red  indicates 

receiver  location,  +  indicates  estimated  location,  and  "o"  indicates  original  location. 


Figure  B.28:  Node  44.  Geometry  during  AEL  event  (top)  and  source 
track  upon  convergence  (bottom). 

Day/h/min:  084/15/40 


AS-00-43 
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